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ABSTRACT 


U 


The supernova acceleration model for cosmic rays of Colgate is 
examined to see if heavy nuclei can indeed survive catastrophic accele- 
ration. If it is assumed that all species of nuclei are collectively 
accelerated by a mechanism such as a plasma wave instability in a 
shock wave, then the only agent remaining to dissociate heavy nuclei 
in Colgate's model is the blue shifted photon flux in the radially 
outward moving shock wave. This would cause a sharp cutoff of the 
heavy component of cosm'.c rays from photodisintegration for cosmic ray 
energies greater than 10 GeV/nucleon. 
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I. INTRODUCTION: 

Colgate and co-workers (I960, 1962, 1963, 1965, 1966, 1967) We 
developed a model for accelerating the mass of the outer mantle of a 
massive supernova (M « 10 Mg) to relativistic energies, thus providing 
a source mechanism for the production of cosmic rays. Their mech- 
anism is based on the formation of a strong shock wave in the super- 
nova process which moves radially outward through the stellar envelope 
with increasing velocity, reading relativistic velocities when a 
mass fraction equal to 10"^ of the total stellar mass remains external 
to the shock. 

The shock wave is formed when the fusion process depletes the 

elements up through Fe in. the core thus causing rapid coding of the 

* 

core. To maintain equilibrium, the inner mantle of the star gravi- 
tationally collapses onto the core, raising the temperature to an 
energy equivalent of kT ft* 50 MeV in a thin outer layer of matter of 
the essentially hard core material. At such high temperatures the 
core material quickly is reduced from a mixture which is predominantly 
Fe 56 , He^, p, n, and e by thermal decomposition to a mixture of p, n, 
and e. The energy is removed from this reaction region by neutrino 
emission by the process p + e -> n +V and redeposited in a thin region 
external to the core by inverse y/?-decay. The redeposition process 

occurs in a thin region external to the core raising its temperature 

$ 

rapidly and thus producing the shock wave. 

In Figure 1 is shown a plot of the shock wave's radial position 
as a function of time, along with the shock wave velocity. The different 
curves are labelled with the mass fraction external to the shock. 
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The temperature versus mass density is plotted in Figure 2, showing 
the neutrino coding of the core and heating of the inner mantle which 
results in the shock wave, as well as the heating of outer mantle by 
the shock wave. 

It is noted that this described mechanism of supernova explosion 
is different from that of Burbidge, Burbidge, Fowler, and Hoyle (1957). 

The latter invoked a thermonuclear explosion as the result of the 
gravitational collapse heating the inner layers of the star. Colgate 
and White (1966) show that the change in gravitational potential energy 
of the imploding core is greater than that of any such thermonuclear 
process. They show, instead, that even should such an explosion occur, 
the rarefaction created by the imploding core material can "swallow" the 
energy of a thermonuclear explosion because the velocity of sound in 
the imploding matter is greater than in the external matter in an 
explosion. 

To obtain the total amount of energy released as cosmic rays by 

a supernova, Colgate and White (1963) consider a 10 star that 

undergoes the supernova process, to be a type II supernova. When the mass 

fraction external to the shock wave is 10”^ , the kinetic energy of an 

accelerated proton is equal to its rest mass. They compute then, to 

a reasonable approximation, that the total energy produced in the form 

of cosmic rays is just 10“^ x 10 M c^ = 1.8 x 10-^ ergs, for a single 

0 

type II supernova. Taking the energy density of cosmic rays in the 
galaxy to be 5 x 10“ ^ erg/cm^, the galactic volume equal to 5 x 10^ cup 
and the cosmic ray escape time as 2 x 10^ yrs, would require one type II 
supernova every 10^ years. Using a higher estimate for the galactic 
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cosmic ray energy density of 10“^ erg/cnr* (Ginzburg and Syrovatskii, 1964) 
would require a higher frequency of type XX supernovae of one every 10^ 
years. These frequencies are commensurate with the currently estimated 
frequency of supernovae in our galaxy, hence the model is energetic 
enough to qualify as the source of cosmic rays. 

The hydrodynamic calculations have not been explicitely extended 
into the relativistic region, but Colgate and Johnson (1960), used 
similarity solutions to compare the analytic one-dimenSional problem 
to the three dimensional case of the supernova. Making the assumption 
that the material behind the shock has uniform energy per rest mass 
they found the proportionality 

yLc CC 

where 4 ~ 7 -I is the kinetic energy per unit rest mass and is the 
mass density of the material. Allowing for the variation of energy 
per rest mass they obtained the lower limit of 


-/ 


cc, //ft 

Withya and T designating the density and temperature, respectively, at 
the shock the mass external to the shock is 


M (< /i ) cc /oT. 

From their hydrodynamic calculations 

T cc/^- 3/ 

so that the number of particles with energy greater than ju is given by 
N (^ > / U ) / V i (^/°) CC yO ^ & 1 

Using the upper and lower limits for kinetic energy as a function of 


density in this last expression they obtain 
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i j / \ -^ < n - 1*${ 

N (>u) oc ^ 

which is in good agreement with the observed cosmic ray integral 
spectral index of n « 1.5. 

The shock conditions leading to the acceleration of the stellar 
matter to cosmic ray energies in Colgate’s model are essentially 
as follows. In the shock wave itself the energy equivalent of the 
temperature is of the order of IQ 2 KeV. The energy density in the 
shock wave region may be written as 

e * RT + aT y /<} , 

Using values for ^ and T from Figure 2 it is easily seen that 



indicating that the energy in the shock wave is almost entirely 
contained in the form of electromagnetic radiation with a black-body 
distribution of temperature T in the shock rest frame. Using the Wien 
displacement law for unit frequency, 

V. „ T = cT/v „ - o.f/e* a K 

one finds that 200 KeV is equivalent to a temperature of = 10^ °K. 

In order to set an uppeir limit on the energy obtainable from such 
a radiation shock wave, Colgate and White (1963) note that the 
condition for the containment of the radiation is that the matter 
external to the shock must be a mean free path for the radiation or 
about 1 gm/cm^. For their stellar model, setting 


J /> c/i” ~ / <jm fan * “ f 




with the scale height h 0 « 6 x 10^ cm, they find the minimum density 


to be 


~<s? /. *9 

^ , - 8 K /O fyn/din . 

fy»i > i 7 ' 


From the previous expressions for the energy spectra and mass density 
relations they find that E max = 3 x 10 7 GeV for the maximum energy 
that is obtainable with the above described shock wave. This is not 
in disagreement with current observations and the accepted view for 
cosmic rays of galactic origin. 

One remaining critical aspect of Colgate's model is its agreement 
with cosmic ray charge spectra. Ginzburg and Syrovatskii (1964) 
reject any catastrophic mechanism, such as Colgate's, for the simultaneous 
rapid acceleration of both protons and heavy nuclei, their reason being 
that fragmentation of the heavy nuclei will occur on the proton 
component and specifically in Colgate's model the heavy nuclei will 
be destroyed through phot<Pdisintegration by the high photon flux in 
the shock wave. To better understand this problem, the microscopic 
aspect of particle acceleration in the shock wave must be considered. . 

The mechanism presented by Colgate (1965) is .acceleration of the 
protons and heavy ions in the electric field set up by the Compton 
scattering of the electrons ahead of the ions. This occurs since the 
radiation in the shock "seen" by the unaccelerated ions is blue shifted 
into the Compton region. Effectively then the shock carries along 



a concentration of electrons which In turn accelerate the ions. 

Colgate (1965) also suggests that plasma ion waves resulting 
from the counter-streaming of ions of different charge- to-mass ratio 
could accelerate the different species of ions collectively. Two 
alternatives for preservation of heavy nuclei in ?. he acceleration 
process are possible: (1) the original nuclei present in the pre- 

supernova stellar envelope are preserved intact through the shock 
acceleration, and (2) the .heavy nuclei are destroyed by fragmentation 
and photodisintegration and then are re- synthesized by the neutron 
capture process 

(A, Z) + N (A . + 1, Z) + y 

(A + 1, Z) + P (A + 1, Z + 1) + N 

in the high temperature neutron rich plasma immediately behind the 
shock. Colgate has determined that rapid cooling prevents the latter 
process from producing elements beyond Z & 30. It is necessary 

then to depend on the first alternative and preserve the charge 
spectrum through the acceleration process. 
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XX. PHQTQD IS INTEGRATION OF HEAVY NUCLEI 

A simple view of the charge separation in the shock is illustrated 
in Figure 3, with the electrons being Compton scattered toward the 
leading edge of the shock resulting in an electric field parallel to 
the shock velocity. 

It has been noted ^y Colgate (X965) that the charge separation 
acceleration of the stellar material will lead to large relative 
velocities between protons and other nuclei whose charge-to-mass ratio 
is of the order of zfz/A - ( /%. . After shock passage this effect will 
cause fragmentation of the heavy nuclei on the proton component with 
a probability of survival of the heavy nuclei of the order of lO"^, 

As mentioned in the previous section, it is desirable then to have 
the acceleration of all components occur at the same rate, so that 
there will be no relative velocity produced between the components. 

Colgate suggests then that any initial Counter-streaming set up between 
the heavy nuclei and protons will produce a two stream plasma instability, 
and this in turn will accelerate all components to the same velocities. 

If indeed the latter /process will occur, then spallation is no longer 
a consideration in the survival of the heavy nuclei in the shock wave 
acceleration model. 

However, there remains the problem of the photodisintegration of 
heavy nuclei in the blue shifted proton flux of the shock wave. The 
experimental threshold for the giant electromagnetic resonance is in 
the range of 5 to 10 MeV for the incident photons. This limit is 
well exceeded in the region where the stellar matter is being accelerated 
to cosmic ray energies . 

A semi-quantitative theory of the nuclear giant resonance is 
discussed in several references (Blatt and Weisskopf, 1952; DeBenedetti, 1964). 
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This resonance, primarily resulting from the electric dipole transition, 
can cause the emission of a nucleon if the nucleus absorbs a photon of 
energy greater than the average binding energy per nucleon. The most 
probable reaction is the absorption of a photon and emission of a 
single neutron, 

11 ^ (I', A/)* M A 

For most heavy nuclei the resonance peak occurs at approximately 20 MeV, 
and the total cross-section for absorption of a photon at this peak 
energy is given by 


(T 

WAX 


^ JL 
- 3 


A 






-27 


Chi 


( 8 ) 


f 


where A is the mass number, 

A schematic representation of a typical giant resonance for an 
intermediate mass nucleus is shown in Figure 4. For the purposes of 
the present discussion it is convenient to fit a Gaussian to this 
curve so that integration over the cross-section with respect ,to 
energy may be facilitated. An expression that gives a reasonably good 
fit is 



f t 

where t'lV is the energy of the incident photon in the frame of reference 
of the nucleus, E m = 22 MeV, and r = 2.55 MeV. Over most of the 



range of nuclear mass E m will vary by only about 10% from the value 
given, hence the above is quite adequate for the purposes discussed 
here . 


The electrons and nuclei see the same flux of photons in the 
shock wave, hence the photon flux that accelerates the nuclei through 
their coupling with the electrons, are also capable of causing photo- 
disintegration of the nuclei. If 5^ denotes the photon flux seen 
by the electrons and nuclei, and CQ is the total cross-section for the 
absorption of al photon by a nucleus, then 


dr,j = Oj dt < 10 > 

is the number of * photodisiutegrations that a nucleus will undergo in 
a time dt, during the acceleration process. 

In principle it is only necessary to integrate (10) over the 
time of passage of the shock wave past a single nucleus to find the 
total number of disintegrations occurring. This, of course, requires a 
detailed knowledge of the shock wave structure. To make such a calculation 
more, or less, model independent it is desirable to make as few 
assumptions as possible. The approach to be taken here assumes first 
that Colgate’s acceleration mechanism can indeed accelerate the stellar 
matter to cosmic ray energies, and secondly that, the shock wave is 
coupled to the nuclei through the Compton scattering of the photons off 
the electrons and these in turn drag along the ions, either by means of 
a charge separation field, or by the creation of a two stream plasma 
instability. 



The force on a single electron caused by the photon pressure is 


given by 



( 11 ) 


where is the Compton cross-section and is the momentum 

transferred to the electron by a single photon. Now if an ion is coupled 
to the electron by the charge separation field only the equation of 
such an ion is given by 



( 12 ) 


On the other hand if the ion is accelerated collectively because of 
the plasma instability, (12) must be rewritten as 



(13) 


where Z and A represent the nuclear charge and mass number respectively 
and is the momentum of the nucleus. 

Eliminating dt between (10) and (13) gives > 


Ay. _ _ JL £l i£L 

dy 'y- ^ <r e Afr 


(14)' 


with % cancelling as well. The'^total number of photodisintegrating 
collisions that a nucleus undergoes is then obtained by integrating (14) 


from the initial momentum of the nucleus in the shock frame, 


Pi 


to 
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the final momentum as it leaves the shock, , or 


7l J =jr f 1 ' £kL ^ 
J A J n <Tc AP, 

Y% 


(15) 




To evaluate the integral in (15), first consider the Doppler 
shift in the photon energy as seen by the electrons and nuclei, 


hy's, COS&O ItV (l + -ItOS £')/(!" S) ( 16 ) 


with primed quantities referring to the rest frame of the electrons 
and nuclei and the unprimed quantities to the shock rest frame. A 
plot of (16) is shown in Figure 5, indicating that for/*? 0 . 9, the 
blue shifted radiation seen by the nuclei is confined to a core of 
&*•< 30°. Since for greater angles than this the radiation is red- 
shifted to very low energies, it is justifiable to assume that the 
radiation seen by the particles may be treated as if it were a plane 
wave for/? >0.9. In this approximation (16) reduces to 


h V > ~ 2/ hv 


( 17 ) 


where y is the Lorentz factor. A further approximation is to treat 
the photon flux as if it were monoenergetic , so that hv represents 
the peak energy of the black-body radiation in the shock frame. 
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For a nucleus of mass AM, the differential momentum may be 
expressed in terms of the Lorentz factor as 

J ft { r i. 0" M c Yciy (is) 


The Compton total cross-section is expressed in general by the 
Klein-Nishina formula. For large y however it is sufficient to use the 
extreme relativistic asymptotic limit of this formula (Jauch and 
Rohrlich, 1955); 



•A-' 








(19) 


with Ojr = 665 mb being the Thomson cross-section, and the dimensionless 

( 

variable is given by 





( 20 ) 


with 




( 21 ) 


The energy of the scattered photon in the Compton process can be 


expressed as 



( 22 ) 


h V — hv J /[l~h (Ity t*~ to'b&Qj 

» 

with denoting the scattering angle. From this one may express 
the net change of momentum as 




- CV>5&0/f / 4 f«i©'))7(23) 


A reasonable approximation is to take the intermediate value when 
cos ~ 0 , so that 


(zrkv/c)(.zdr)/(HZc<'r) ( 24 ) 

Substituting from (9), (18). (19), and (24) info (14) and using 
(17) one obtains for the total number of photodisintegrations undergone 
by a nucleus in the acceleration process 


X 


i A i/a 

CT- Uv 


(14 2 < ciT)eKf> t~ i 

(r C 4 ir 4 •/ 1. ) 


(25) 


where and y^ represent the initial and final Lorentz factors of the 
nucleus with respect to the shock wave, respectively. When the nucleus 
enters the shock wave it has a velocity equal in magnitude to the shock 
wave velocity, i.e., 7 ^ = y Q . There is a singularity in the integral 

of (25) at 7 = 1 which reflects the fact that an infinite time would 
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be required for the particles to be completely stopped in the shock 
frame. In other words, the particles leave the shock with some residual 
velocity corresponding to y > 1, in the shock frame, and the particles 
immediately after shock wave passage will have a velocity in the rest 
frame of the star of 7 < y Q , 

Above if (12) had been substituted into (10), one would have 
obtained ultimately for (14), 



(26), 


for the case of the charge separation field acceleration. It is 
possible then to express the result for the two cases as 



‘t'd, f 


with 



C 




I AU X 

6 ii V 


(27), 


(28), 


(2.9), 





-16- 


and 
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z f 


(30). 


f( 7 ) is just the integrand in (25). 

In Figure 6 , f( 7 ) is plotted as a function of 7 , showing the 
resonance peak at 7 ® 55 and the singularity at 7 = 1. It is 

easily seen that the lower limit of integration 7 ^ may be set over a 
rather wide range of values such that yf 1 without making any 
substantial difference in the result of the integral. Therefore 7 ^ 
can safely be set equal to 2 without changing the final result appreciably. 
This means only that the final velocity of a proton in either case, 
or a heavy nucleus in the plasma instability case will be equal to 
7 g - 1 in the rest frame of the star where 7 g i’s again the shock 

velocity. ' 

The integral occurring in (27) and (28), is easily performed by 
numerical integration. From the resulting number of collisions in 
each case the probability of survival of the heavy nuclei is simply 

P(survival) - e ^ (31) 


The calculations have been done for iron, with (5^ = 665 mb, 

Mc^ - .938 GeV, and hv = 


.2 MeV. The results for both cases are 
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tabulated in Table X for a range of values of the shock wave velocity 

7 g , These same results are plotted in Figure 6 . 

Xt is readily seen from Figure 6 , that there is little difference 

between the two acceleration mechanisms, but that both show a steep 

cutoff of heavy nuclei at approximately y Q * 35. Along the top of 

the plot are scales showing the corresponding momentum per nucleon 

to which such particles would be accelerated by the shock wave passage. 

For the case of field acceleration this corresponds to about 15 GeV/ 

nucleon, and for plasma wave acceleration to 32 GeV/nucleon. 

Colgate and White (1966) noted that in the non-relativis tic case 

the heated material behind the shock would expand adiabatically after 

shock passage and that its bulk velocity would double in the process, 

corresponding to an increase in kinetic energy by a factor of four. 

Colgate (1967), however, shows that in the relativistic case the final 

energy of the stellar matter after expansion would be proportional to 
o 

7 . Taking this latter result into account, the above cutoff energies 

s 

for the complete destruction of heavy nuclei would occur at 225 GeV/ 
nucleon for the case of the field acceleration process, and 1024 GeV/ 
nucleon for the plasma wave acceleration process, 
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XXX. DISCUSSION 

On the basis of the approximate treatment described in the 
previous section, Colgate's model permits acceleration of heavy 
nuclei up to quite large final kinetic energies. Since Colgate (1968a) 
has shown that it is not likely that the heavy nuclei can be re- 
synthesized in the hot material behind the shock, if spallation 
were to occur, it seems, reasonable to rule out the possibility of the 
field acceleration process since it would surely produce spallation. 

One is left then with a plasma wave process, or some other mechanism, 
which will collectively accelerate all nuclei to the same velocities, 
Such a process seems to allow acceleration of the supernova material 
up to kinetic energies of the order of 10^ GeV/nucleon before photo- 
disintegration on the photon field of the shock wave eliminates the 
heavy nuclei component. 

At the present time experimental observations extend only as far 
as about 50 GeV/nucleon (von Rosenvinge and Webber, 1968) and these 
results show no reduction in heavies. The conclusion that is reached 
then is that so far Colgate's supernova acceleration is not in any 
appreciable disagreement with current cosmic ray observations. 

It should be noted that the version of the model presented 
above is perhaps an extreme case. Since most supernovae are 
probably of mass < 10 M , this may require that some small modi- 
fications be made to the parameters. Further, according to Colgate 
(1968b) there may be a good possibility for resynthesis of the heavy 
nuclei following shock passage by means of successive rapid neutron 



captures 

neutrons 
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i- j 
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followed by p, n charge exchange reactions between bound 
and free protons. 
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Figure Captions 


Figure 1: 


Figure 2: 


Figure 3: 


Figure 4: 


Figure 5: 


Figure 6: 


10 Mq supernova radius versus time, 
(Colgate and White, 1966). 

10 M supernova temperature versus 
density with neutrino deposition, 

(Colgate and White, 1966). 

Schematic illustration of the charge 
separation of ions and electrons in the 
shock region of a supernova, resulting 
in an accelerating field 'fi. 

Schematic illustration of qualitative 
theoretical expectation of the cross- 
section for formation of an excited 
nucleus by gamma-ray absorption for an 
intermediate nucleus, (Blatt and 
Weisskopf , 1952) . 

The angular dependence of Doppler 

shifted radiation for various values 

of /s \ 

/ 

The function f(y) versus y s and the 
probability of survival of an iron 
nucleus for charge separation field 
acceleration and plasma wave acceleration. 
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